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The set of differential equations obeyed by the redshift in the general /3' 7^ Szekeres spacetimes 
is derived. Transversal components of the ray's momentum have to be taken into account, which 
leads to a set of 3 coupled differential equations. It is shown that in a general Szekeres model, and in 
a general Lemaitre - Tolman (L-T) model, generic light rays do not have repeatable paths (RLPs): 
two rays sent from the same source at different times to the same observer pass through different 
sequences of intermediate matter particles. The only spacetimes in the Szekeres class in which all 
rays are RLPs are the Friedmann models. Among the proper Szekeres models, RLPs exist only in 
the axially symmetric subcases, and in each one the RLPs are the null geodesies that intersect each 
t — constant space on the symmetry axis. In the special models with a 3-dimensional symmetry 
group (L-T among them), the only RLPs are radial geodesies. This shows that RLPs are very 
special and in the real Universe should not exist. We present several numerical examples which 
suggest that the rate of change of positions of objects in the sky, for the studied configuration, is 
10~® — 10~^ arc sec per year. With the current accuracy of direction measurement, this drift would 
become observable after approx. 10 years of monitoring. More precise future observations will be 
able, in principle, to detect this effect, but there are basic problems with determining the reference 
direction that does not change. 



THE MOTIVATION 



The quasi-spherical Szekeres solutions have recently 
begun to be taken seriously as cosmological models [l[ - 
[31 • For this application, one has to know the equations 
obeyed by the redshift. The corresponding equation for 
radial null geodesies in the Lemaitre - Tolman (L-T) 
model [1, Q was derived long ago by Bondi [l^ , see also 
Ref . [ll| . The generalisation to the Szekeres geometry is 
nontrivial because in general there are no radial geodesies 
in the latter [ij, [l3 ■ Consequently, the transversal com- 
ponents of the ray's momentum necessarily have to be 
taken into account, and a set of 3 coupled differential 
equations is obtained. These equations can then be ap- 
pHed to nonradial geodesies in the L-T model. 

The purpose of this paper is to derive the redshift 
propagation equations in a general^ Szekeres model of 



the /?' ^ family [ll[, so that they can be numerically 
solved and applied in various situations. 

In Sec. [ni the Szekeres models are introduced. In Sec. 
mil it is pointed out that the Bondi redshift equation for 
radial null geodesies in the L-T model is in fact an ap- 
proximation, the small parameter being the period of the 
electromagnetic wave. The same is true for the equations 



derived here. In Sec. llVi the general equations of null 
geodesies in Szekeres models are presented. In Sec. El 
the set of redshift equations for the Szekeres models is de- 
rived. In Sec. IVH conditions are discussed under which 
light rays between a given source and a given observer 
proceed through always the same intermediate matter 
particles; such rays are termed "repeatable light paths" , 
RLPs. In Sec. IVIH the equations of Sees. Ivl and IVTl are 
applied to general null geodesies in the L-T model and in 
the associated plane- and hyperbolically symmetric mod- 
els. It is shown there that in these models the only RLPs 
are the radial null geodesies. Sec. llXl is a brief summary 
of the results. 



II. THE SZEKERES SOLUTIONS 



The Szekeres solutions [Tj, 1T5| follow when the metric 



ds' = dt' 



,2a(t,r,.,y)^^2 _ g2/3(t,r,.,,;) (^^^2 ^ ^^2^ 
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The generalisation to 
cover the quasi-plane and quasi-hyperbolic cases is immediate, 
so it would not make sense to leave it out. 



(2.1) 

is substituted in the Einstein equations with a dust 
source, assuming that the coordinates of ()2.1|) are co- 
moving, so that the velocity field is u^ — 6^q (with 
{x" , x^ , x'^ , x^) = {t,r,x,y)). 

There are two families of Szekeres solutions, depend- 
ing on whether (5,^ = or P,^ ^ 0. The first family is a 
simultaneous generalisation of the Friedmann and Kan- 
towski - Sachs [16] models. Since so far it has found no 
useful application in astrophysical cosmology, we shall 
not discuss it here (see Ref. [li|). After the Einstein 
equations are solved, the metric functions in the second 



family become 

Qp = $(t,r)e'^('''^'^\ 

e« = /i(r)$(i,r)/3,, = /i(r) ($,,+$!.,, ), (2.2) 

e""^ = A{r){x^+y'^)+2Bi{r)x + 2B2{r)y + C{r), 

where the function $(t, r) is a solution of the equation 



*,t = -k{r) 



2M(r) 1 
$ 3 



A$2; (2.3) 



while h{r), k{r), M{r), A{r), Bi{r), B2{r) and C{r) are 
arbitrary functions obeying 



dcf 



7(r) = 4 (AC - Bi^ - 52^) = l//i2(^) + k{T). (2.4) 
The mass density in energy units is 
(2Me3'') ,^ 



K/9 : 



32^^ (e^) ,. ' 



K = SttG/c'^. (2.5) 



Whenever (e^) ,r = and {2Mc^'') ,r^ 0, a shell cross- 
ing singularity occurs. It is similar to the shell crossing 
singularity in the L-T models, but with a difference. In 
a quasi-spherical model a shell crossing may occur along 
a circle, or, in exceptional cases, at a single point, and 
not at a whole surface of constant t and r, as was the 
case in the L-T models. 

As in the L-T model, the bang time function follows 
from (lOll: 



d$ 



-k + 2M/$ 



iA$2 



= t-tBir), 



(2.6) 



The solutions of the above equation for A 7^ involve 
elliptic functions and were first studied by Barrow and 
Stein-Schabes [17j]. 

As seen from (|2.ip and (|2.2p . the Szekeres models are 
covariant with the transformations r — f{r'), where /(r') 
is an arbitrary function. 

The Szekeres metric has in general no symmetry, 
but acquires a 3-dimensional symmetry group with 2- 
dimensional orbits when A, Bi, B2 and C are all constant 
(that is, when iy,r ~ 0). 

The sign of g{r) determines the geometry of the 2- 
surfaces of constant t and r (and the symmetry of the 
constant A, Bi, B2 and C limit). The geometry of 
these surfaces is spherical, planar or hyperbolic (pseudo- 
spherical) when g>0,g = 0or(7<0, respectively. 
With A, Bi, B2 and C being functions of r, the surfaces 
r = const within a single space t = const may have dif- 
ferent geometries, i.e. they can be spheres in one part of 
the space and the surfaces of constant negative curvature 
elsewhere, the curvature being zero at the boundary. 

The sign of k{r) determines the type of evolution; with 
fc > = A the model expands away from an initial sin- 
gularity and then recoUapses to a final singularity, with 



fc < = A the model is either ever-expanding or ever- 
collapsing, depending on the initial conditions; fc = is 
the intermediate case corresponding to the 'flat' Fried- 
mann model (fc = can also occur on a 3-surface as 
the boundary between a region with fc > and another 
one with fc < 0). The sign of k{r) influences the sign 
of g{r). Since 1/h^ in (|2.4p must be non-negative, ^ we 
have the following: With 5 > (spherical geometry), all 
three types of evolution are allowed, with 5 = (plane 
geometry) , k must be non-positive (only parabolic or hy- 
perbolic evolutions are allowed), and with g < (hyper- 
bolic geometry), k must be strictly negative, so only the 
hyperbolic evolution is allowed. 

The Friedmann limit follows when $(i, r) — ^i{r)S{t). 
No further specialization of the Szekeres functions is 
needed; the limiting Friedmann model is represented in 
the little-known Goode - Wainwright [T^ coordinates, 
see also Ref. [l9{ . 

The Szekeres models are subdivided according to the 
sign oi g{r) into the quasi-spherical (with g > 0), quasi- 
plane {g = 0) and quasi- hyperbolic ones {g <0). Despite 
suggestions to the contrary made in the literature, the 
geometry of the latter two classes has not been investi- 
gated at all and is not really understood; see Refs. [20] 
and [21] for recent work on their interpretation. Only the 
quasi-spherical model has been rather well investigated, 
and found useful application in astrophysical cosmology. 
However, including g < in the redshift equations causes 
no complication, so we consider here an arbitrary g. 

The quasi-spherical model may be imagined as a gen- 
eralisation of the L-T model in which the spheres of con- 
stant mass are made non-concentric. The functions A{r), 
Bi{r) and B2{r) determine how the centre of a sphere 
changes its position in a space t — const when the ra- 
dius of the sphere is increased or decreased [22|. Still, 
this is a rather simple geometry because all the arbitrary 
functions depend on one variable, r. 

It is often convenient to reparametrise the Szekeres 
metric as follows [2j] . Even if j4 = initially, a transfor- 
mation of the (x, j/)-coordinates can restore A ^ , so we 
may assume A ^ Q with no loss of generality [ll[ . Then 
let g ^ 0. Writing A = ^\/{2S), B^ = -^\P/{2Sl 

B2 = -/MQ/(25), e = g/\gl k = \g\k and $ = /M*, 
we can represent the metric (|2.2p as^ 



dcf (X-P)^ 

25 



{y-Q? eS 
25 2 ' 



(2.7) 



^ l/h(r) can be zero at isolated points - it is then either a coordi- 
nate singularity or a neck or belly — but not on open intervals. 

^ The tildes were dropped in I I2.7I I and in all further text. The <& 
in 1 I2.7I I is in fact <E> and the k{r) is k{r). The redefinitions imply, 
via lHHl, C = vial [{P^ + Q^) /S + eS] /2, h^ = l/[\g\{e - k)] 



and M ■■ 



\g\ M. The M used from now on is in fact A/. 



ds' 



dt' 



(^,,-$g,,/g)^ 

e — k{r) 



dr^ 



$2 



(da;2 + dy2 



where, so far, e — ±1 (+1 for the quasi-spherical and —1 
for the quasi-hyperbohc model). When 3 = 0, the tran- 
sition from (123) to (1221) is A = 1/(25'), Bi = -P/{2S), 
B2 = -Q/{2S) and $ is unchanged."* Then ([2Jl) applies 
with e = 0, and the resulting model is quasi-plane. 

The parametrisation introduced above makes several 
formulae simpler, mainly because the constraint ()2.4p is 
identically fulfilled in it. However, this parametrisation 
obscures the fact, evident in (I2.ip - (I2.4p . that the same 
Szekeres model may be quasi-spherical in one part of 
the spacetime, and quasi-hyperbolic elsewhere, with the 
boundary between these two regions being quasi-plane; 
see an explicit simple example in Ref. [20]. In most of 
the literature published so far, these models have been 
considered separately, but this was either for purposes of 
systematic research, or with a specific application in view 
that fixed the sign of g{r). 

Equation (j2.3p . is formally identical to the Friedmann 
equation, but with k and M depending on r, so each 
surface r = const evolves independently of the others. 
The solutions $(i, r) are the same as the corresponding 
L-T solutions, and are unaffected by the dependence of 
the Szekeres metric on the (x, y) coordinates. 

As defined by (|2.2p - (|2.3p , the Szekeres models contain 
8 functions of r, of which only 7 are arbitrary because 
of (|2.4p . The parametrisation of (|2.7p turns g{r) to a 
constant parameter e, thus reducing the number to 6. 
By a choice of r (still arbitrary up to now), we can fix 
one more function (for example, by defining r' = M{r)). 
Thus, the number of arbitrary functions that correspond 
to physical degrees of freedom is 5. 

In the following, we will represent the Szekeres solu- 
tions with /3,r 7^ in the parametisation introduced in 
(|2.7p . The formula for mass density in these variables is 



up 



2{M,r-3M£,r/£) 

The shear tensor is 

1 /$,fr-*,t$,r/$ 



(2.8) 



/3 



3 V $,r-*f,r/f 

and the scalar of expansion is 

U = U -a = —T 1 -t-t: 77;— 



diag(0,2,-l,-l), (2.9) 



(2.10) 



III. REMARKS ON THE BONDI REDSHIFT 
EQUATION IN THE L-T MODEL 

The L-T model is a special case of the quasi-spherical 
Szekeres models that follows from (12.71) when e = -1-1 



* The implied changes in C and h are then C = (P^ + Q'^)/(2S), 
h? = — 1/fc; k and M remain unchanged. 



and the functions P, Q, S are all constant. With a dif- 
ferent representation of the coordinates on a sphere, the 
resulting metric is: 



ds' = dt' 



R,, 



-dr^ ~ R^{t, r) (di?^ + sin^ Mip^) , 

(3.1) 



1 + 2E{r) 

and the equation of an incoming radial null geodesic is 

R,r{t,r) 



dt 
dr 



VTT2E' 



(3.2) 



Bondi's derivation [10] of the redshift equation for this 
geodesic is as follows. Take a light signal obeying p.2p . 
the equation of its trajectory (the solution of p.2p ) is 



t = r(r) 



(3.3) 



Take a second light signal, emitted from the same radial 
coordinate r, but later (as measured by the time coordi- 
nate t) by r. The equation of its trajectory is: 



t^T{r) + T{r), 
where (T + t) obeys, from l\'S.2\\ : 

R,r{T{r) + T{r),r) 



dr 

dr 



dr 
dr 



(3.4) 



(3.5) 



v/1 + 2^(0 
From the Taylor formula we have: 

i?„(r(r)+T(r),r) = i?„ (r(r), r) + r(r)i?,,, (r(r), r) 



0{t^ 



(3.6) 



where the last term has the property Oir^^r^/T — > 0. 

Now, assuming that t is sm.all, we neglect the last term 
in (|3.6p and obtain from p.Sp , taking into account p.2p : 



dr ^ R,tr{T{r),r) 

= —T(r) — , -. 

dr Vl + 2^W 



(3.7) 



If T is the period of an electromagnetic wave, then by 
definition: 



T(robs) 

T(rcm) 



1 + ^(''Gm), 



(3.8) 



where the subscripts 'obs' and 'em' refer to the points 
of observation and emission, respectively, and z is the 
redshift. From (|3.8p . keeping the observer at a fixed 
position and letting rem vary, we obtain (dr/dr)/T = 
-(dz/dr)/(l + z), and so in (jXTl) : 



1 



1 



dz 
z dr 



R,tr(T{r),r) 
v/1 + 2i?(r) ■ 



(3.9) 



This is Bondi's radial redshift equation jl3]- It does not 
describe the redshift propagation exactly. Neglecting the 
last term in (13.61) we have changed the exact equation 
into one that only approximates the actual variation of 



T along the ray. The approximation is better the smaller 
the value of t. Considering that r is the period of an 
electromagnetic wave, and taking into account the pe- 
riod range of relevance in observational astronomy (from 
gamma rays up to radio waves, the longest observed of 
which have the wavelength of the order of 15 m, thus 
the period of about 5 x 10'^ s), we see that, compared 
to cosmological time-scales, the periods are short indeed 
and the approximation is not bad. Moreover, as seen 
from (|3.8p . by following the rays back from the observa- 
tion event into the past, we encounter ever smaller values 
of T, so the approximation gets progressively better with 
increasing redshift (or, rather, gets progressively worse 
as the ray approaches us). Still, it is conceptually impor- 
tant to remember that ()3.9p involves an approximation 
(this approximation is equivalent to the geometric optics 
approximation [11, 24] that leads to the commonly used 
expression for the redshift 1 + z — {kaU°')cm/ika'u°')ohs)- 
We shall apply the same approach to the redshift equa- 
tions in the Szekeres models in Sec. El 



IV. EQUATIONS OF GENERAL NULL 
GEODESICS IN A SZEKERES SPACETIME 

For reference, the equations of general null geodesies in 
a Szekeres model are copied from Ref. [l2| in Appendix 
[K\ ^ They are written there in terms of an afHne param- 
eter s. For our present purpose it is more convenient to 
use the coordinate r as an independent parameter (which 
is non-afhne). 

This is allowed, but with some caution. It is easily 
seen from (jAip - (|A4p in Appendix |X] that a geodesic 
on which dr/ds — over some open range of s has 
dx/ds ~ dy/ds = in that range, and so is timelike. 
However, (jAip - (jA4P do not guarantee that dr/ds ^ 
at all points; isolated points at which dr/ds — can 
exist. Examples that explain how this can happen are 
the non-radial geodesies in an L-T model, considered in 
Sec. IVIIII Thus, r can be used as a parameter on null 
geodesies only on such segments where ds/dr > or 
ds/dr < throughout. 

Several sub-expressions in the equations of Appendix 
El are multiply repeated, therefore we introduce the fol- 
lowing abbreviations: 



dof 



dof 



$,t^-$,t£,r/£ 



$01, 



<i>£,rr/£ =^ $11, 



(4.1) 
(4.2) 
(4.3) 



^ It is shown in Ref. [131 (see also Ref. |13|| ) that in general there 
exists no analogue of a radial null geodesic. Radial geodesies 
exist only when the Szekeres model is axially symmetric; then 
their intersections with every space of constant time coordinate 
lie on the axis of symmetry. 



^■jr^^x ^^:xr — -l^12i 
dcf 



(4.4) 

S,r£,y^££,yr = Ei^. (4-5) 

In addition, the following replacement will appear useful: 



dxWfdy_Y drf^_ 



dr y \dr 
We have, for any coordinate 

2^2„ 



ds^ \ds J dr^ ds^ dr 
Then, from (IA2P we have: 



(4.6) 



(4.7) 



ds2 \dsj { $1 dr \^i £ 2e-k 

_ 2— — — -2— — — -H— ^^E 
£2 $^ dr £2 (j)^ dr £2 (j)^ 

'^^Uit,r,x,y)f^V 



(4.8) 



Consequently, (KH . (K3^ and (|X4l) become, using (|i?71) : 
d2t 'I'l'&oi ^^,t dt 

d^^ + T^ + ^^ + ^d^-O' (4-^) 



a2 

dr2 



d^x „<i>,t dt dx 1 $1 






$ dr dr ^ e — k 



Ei-2 



2$i da; _ f^ /da;\ £,y dx dy 

$ dr £ \dr J f dr dr 



(4.10) 



dr2 



^^^didy 
$ dr dr 



1 $1 



-E 



13 



2^dy £yy_ /dx 
$ dr £^ \dr 



-2 



£,x dxdy 
£ dr dr 



£ \ dr J dr 



(4.11) 



V. THE REDSHIFT EQUATIONS IN THE 
SZEKERES MODELS 

Consider, in the Szekeres metric (|2.7p . two light sig- 
nals, the second one following the first one after a short 
time-interval r, both emitted by the same source and ar- 
riving at the same observer of coordinates {r,x,y). The 
equation of the trajectory of the first signal is 

it,x,y)^{T{r),X{r),Y{r)), (5.1) 

the corresponding equation for the second signal is 

(i, X, y) = {T{r) + T{r),X{r) + C(r), F (r) + ^(r)). (5.2) 



This means that while the first ray intersects the hy- 
persurface of a given constant value of the r-coordinate 
at the point {t,x,y) — {T,X,Y), the second ray in- 
tersects the same hypersurface at the point (i, x, y) = 
{T+T,X + (, Y + ip). Thus, in general, those two rays will 
not intersect the same succession of intermediate mat- 
ter worldlines on the way. Note that the coordinates 
we use throughout the paper are comoving, so both the 
source of light and the observer keep their spatial coor- 
dinates unchanged throughout history. Given this, and 
given that we consider a pair of rays emitted by the 
same source and received by the same observer, we have 
{C,tp) = (0,0) at the point of emission and at the point 
of reception. However, we have to allow that the second 
ray was emitted in a different direction than the first 
one, and is received from a different direction by the ob- 
server.^ The directions of the two rays will be determined 
by {dx / dr , dy / dr) and {dx/dr + S,(r),dy/dr + r]{r)), re- 
spectively, where ^ = dC/dr, 77 — dtjj/dr. We will assume 
that {dT/drX,'fp,^,il) are small of the same order as r, 
so we will neglect all terms nonlinear in any of them and 
terms involving their products. 

Since ^ = V-" = at the observer, these quantities are 
not in fact observable. However, they have to be numer- 
ically monitored along the ray because, as will be seen 
below, they enter the equation for r, which is connected 
to the redshift by ([3^. 

In writing out the equations of propagation of red- 
shift, we will introduce the symbol A. It will de- 
note the difference between the relevant expression taken 
at {t + T,r,x + (^,y + tp) and at {t,r,x,y), linearized 
in (r, C, V'); for example ^Qi{t + T,r,x + (,y + ijj) — 

$oi(t,r,a;,y) ='a$oi+0(t2,tC,t^,C',-.-)- We have: ^ 



A$ = $,tT, A($,t) =$,ttr, 

di dr a'^^ c a "^^^ 

dr dr ' dr ' dr ' 

A£ = £,xC + £,yip, 

AS,, = C/S, A£,y^^/S 



£2 ^ ' £2 



'^,tEi2, , <^,tE 



£2 



c + 



13 



52 



(5.3) 

(5.4) 
(5.5) 

(5.6) 



A$ii = {<^,trr-'^t£rrr/£)T+ — {£,rr£,x-££,rrx)C 



£2 



£2 



~r ^2 V^ ''"^ ^ 'y ^^ trry ) W • 



(5.7) 



^ This means that in a general inhomogeneous and anisotropic 
Universe the observed objects should drift across the sky. See a 
brief quantitative discussion of this effect in Sec. IVIIII 

"^ A quick way to calculate 115.31 1 — I I5.16I I is to take the differ- 
ential of the corresponding quantity at constant r and replace 
(dt, da;, Ay, d(da;/dr), d(dy/dr)) by (t, C, V", ?, v)- 



In the next two equations account is taken of the fact 
tnat c- ,xy -— ^■ 

^-L^12 — \^ jr ^ txx ^^ ^rxx ) ^ ~r \^ ^ry ^ jx ^ :y ^ jrx ) V^^ 

(5.8) 

ZA£/i3 y^ ,rx ^ ty ^ -iX ^ iry } s ~r yO ^r ^ tyy ^^ ^ ^"^yy ) 'r^' 

(5.9) 



A ^ „ da; , „ dw 
AS = 2— e + 2/77, 
dr dr 



(5.10) 



AU = 2{ -^£21 + ^01^^ A ^ _ 2— — 



$1 $^2 J ^j. (j)^ (J^ 

A$ii $iiA$i A£,^ £,rAf 



$1 $1^ 



£ 52 



£^<t>i £3$i £2$^ 

<^Ei2A^i\ dx_^<^Ei2^ 



^,tEi3T J>A£Ei3 $Agi3 

£2$^ £3$^ ^2$^ 

^Ei3A^i\ dy ^^Ei^ri 



£2$^2 J J^ £2$^ 

(£-fc)$S /'$,tT A£ A<I>i ^ AS 



£2$^ y $ £ $1 S 



(5.11) 



Applying the A-operation to (|4.9p - (|4.1ip we obtain: 
d2T $oiA$i +$iA$oi {<i>,t'^ + '^'^,tt)'^T 



dr2 



e — k 



£2 



£'^ c^ dr dr 

(5.12) 



^ + 2^^-^^— — r + 2^ — — 
dr2 \ ^ ^'^ J dr dr $ dr dr 

$,t dt A^iE^ia ^,t^iEi2T ^lAE^ia 
$ dr^ (e-fc)$ (e-/c)$2 (£-fc)$ 
/A<l>i <i>i<l>,t r\ dx f,*i(. 

dxV /^ _ £,xA£ \ _ 2£,, g dx 
dry V-Sf '^^ / £^ dr 

jdx dy / V _ ^^yA^X _ £,y (dy dx 
'drdr \S£ £2 ) ^ c 1 .„4+,„^ 

2 



i5 V dr dr 



dy\ Y^ _ g,.Ag \ 2g,^ 77 dy 
dr) \S£ £^ J £ dr 

AU^ + U^ = 0, 
dr 



ft + o/^,» $,t'\ dtdy_ ^ ^<i>,tdydT 
dr2 V $ $2 / dr dr <i> dr dr 



(5.13) 



<i> dr 



'A$i $i$,tT\ dy $1 
2'^ i^)d^ + 2^^ 



A$iJ;13 $,t$i£:i3T $iAj;i3 
(e-fc)$ (e-fc)$2 (£-fc)$ 



' dr/ IS-f 



g,yAg 

£2 



- 2- 



dx dy /^ C ^!i: ^'^ 

£,yA£ 



drdr V-S^ 



£2 



2£,y£^ dx 
£ d7 

2g,y T? dy 

f dr 



AC/^ + C/r; = 0, 
dr 



da; 
dr 



(5.14) 



In addition, we have the first integral of the geodesic 
equations (|i^ - (|i?TT|) : 



dtV ^ ($i)^ 
dr / £ — /j 



$2 



dx 
dr 



(5.15) 



Applying the A-operation to this we get 
$iA$i 



drdt 
dr dr 



dx 
dr 



e — k 

$^ fdx dy 
f 2 \^ dr dr 

Note: dt/dr < for an incoming ray. 



(5.16) 



VI. REPEATABLE LIGHT PATHS 



As attested by (|5.12p - (|5.14p . in a generic Szekeres 
model two light rays connecting a given source to a given 
observer at different instants of emission do not proceed 
through the same succession of intermediate matter par- 
ticles. We will now investigate under what conditions 
this intermediate succession is the same. This property 
will be called repeatable light paths (RLP). 

For a RLP we have 



C=v=^=^=o 



(6.1) 



all along the ray. Then (|5.12p decouples from (|5.13p - 
(|5.14p and just determines r (and, with it, the redshift), 
if the null geodesic equations are solved first. Equations 
(|5?T3l) - l[5A^ become then: 



{e-k)^ 



*, 



dtdy 



<i>2 J dr dr 

^ $,t$l£'l3T 

(e - fc)$2 
$i£^\ dy 
$2 J dr 



^^^dydr 
$ dr dr 



dr 



0. (6.3) 



These equations can be understood in 2 ways: 

1. As equations defining special Szekeres spacetimes 
in which all null geodesies are RLPs. 

2. As equations defining special null geodesies which 
are RLPs in subcases of the Szekeres spacetimes. 

In the first interpretation, (|6.2I) - (|6.3I) should be iden- 
tities in the components of dx"/dr. They are polynomi- 
als of degree 3 in these components, and since (di/dr)2 
does not appear in them, the constraint (I5.15P plays no 
role - all powers of dx"/dr that do appear are indepen- 
dent. Equating to zero the coefficient of (dx/dr)^ in (|6.2p 
(which arises inside AU, within E), and taking into ac- 
count that A£ — AS — when (|6.ip holds, we get 



(6.4) 



*='$,(. -$,*$,./$ = 0. 



The integral of this is $ = S{t)f{r), where S and / 
are arbitrary functions. It is seen from p.9p that this 
means zero shear, i.e. the Friedmann limit. With 
fulfilled, (|6.2I) and (16.31) become identities, and 
(14. lip reduce to the equations of general null geodesies in 
a Friedmann spacetime.* With the observer placed at the 
origin, the geodesies become radial, dx/dr = dy/dr = 
0, and then (|5.16p becomes equivalent to the ordinary 
Robertson - Walker redshift formula, 1+z = S{to)/S{te)- 
To verify this, some calculations are needed, in which Ref. 
|11| may prove helpful. 

Thus, we have proven the following: 

Corollary 1: 

The only spacetimes in the Szekeres family in which all 
null geodesies have repeatable paths are the Friedmann 
models.^ 

In the second interpretation of (|6.2p - (|6.3p . we con- 
sider 2 cases: 

A. The general case: dx/dr 7^ 7^ dy/dr every- 
v^rhere. 

Then we multiply (j6.2p by dy/dr, (j6.3p by dx/dr and 
subtract the results. Disregarding the familiar case ^I* = 
we get 



dy dx 

-c/12-j -C/13-— — u. 

dr dr 



(6.5) 



(j) $2 y dr dr 

A$ii?i2 ^,t^lEi2T 



(£-fc)$ 



(e - fc)$2 

$l$,tT 



$,( dxdr 
$ dr dr 



$2 



^ + Af/^ = 0, (6.2) 
dr dr 



* We recall, however, that the Friedmann limit is represented in 
the Goode — Wainwright [lS| coordinates (see the remark in para 
4 after 1 12. 6I 1). Consequently, all equations representing the Fried- 
mann model will look unfamiliar. 

® This is one more piece of evidence of how exceptional the Robert- 
son - Walker class of models is. 



This, together with dO]), ((5?T2|) and (gj]) - (|4TT1) defines 
a certain subcase of the Szekeres model and a class of 
curves in it. Since both the subcase and the class will 
turn out to be empty, but the calculations proving it are 
rather elaborate, we present them in Appendix [Bj 

B. The special cases: dx/dr = or dy/dr = 0. 

These two cases are equivalent under the coordinate 
transformation {x,y) — (j/',a;'), so we consider only the 
first one. Again disregarding ^ = 0, we get from (|6.2p 
Ei2 — 0. Then, (I4.10[) implies two possibilities: 

Ba) £„ = 0. 

This is possible only when P is constant, and then the 
geodesic lies in the subspace x = P. Equations (I4.11|) and 
(|01) still have to be obeyed, while (l^lQ]) and dOI are 
fulfilled identically. The simple coordinate transforma- 
tion X = x' + P has then the same effect as if P = and 
X = along the geodesic. We show in Appendix |E] that 
in this case, apart from the axially symmetric subcase 
mentioned below, RLPs may exist only when the Szek- 
eres metric has a 3-dimensional symmetry group. Such 
spacetimes are considered in Sec. IVIIl 

Bb) dy/dr = 0. 

The case da;/dr = dy/dr = 0, £ = +1 was investigated 
in detail in Ref. [12| . It turned out that this can happen 
only when the Szekeres spacetime is axially symmetric, 
and then along only one sub-family of null geodesies ~ 
those that intersect each t = constant space on the sym- 
metry axis. We show in Appendix [F] that this result 
applies also with e < 0, and that other RLPs may exist 
only with higher symmetries. 



(|4TT|) and ((0|l are fulfilled identically. However, ((O 
is not an identity and reduces to: 



da; 
dr 



*,tt *, 






dt „$,tdr 
dr <3> dr 



2 - 



$,* 



*,t 



<I> 



$,r 



dt 



~2— ; T H T 



$,^ dr 



$, 



(e — k)T f dx 
dr 



£2 



dof da: 

= :5- X = 0. 

dr 






<J> 2 



(7.1) 



One solution of this is dx/dr — 0, which together with 
dy/dr — defines a radial null geodesic. Then, (|4.10p ~ 
(|4ll|) are fulfilled identically, while ((5T5)) - ((5lB . to- 
gether with (|4.ip - (|4.2p and (|5.5p reproduce the Bondi 
equation p.7p when e — +1. So, we found that in the 
G3/S2 models all radial null geodesies are RLPs.^° 

There would exist other RLPs in these models if x in 
(17. ip were zero along any null geodesic - possibly in some 
subcases of the models. It is shown in Appendix [Gj that 
this does not happen, so the radial null geodesies are the 
only RLPs in these models. 



VIII. NUMERICAL EXAMPLES OF NON-RLPS 
IN THE L-T MODEL 



VII. RLPS IN THE G3/S2 MODELS 

The symbol G3/S2 denotes such models that have 3- 
dimensional symmetry groups acting on 2-dimensional 
orbits [ll|. They result from the general /?' 7^ Szekeres 
family when the functions (P, Q, S) are all constant. The 
symmetry of the model is then spherical when e — +1 
(this is the L-T model), pseudospherical (also called hy- 
perbolic) when e = — 1 and plane when s — 0. 

Using the G3 symmetry, the origin of the (a;, y) coor- 
dinates at (a;, y) = (P, Q) can be moved to any location 
on the S2 surfaces. So let us consider the ^'2 on which 
the first light ray is emitted, and let us choose the origin 
of (a:, y) at the position of the emitter. Thus, in (|4.9p 
- (|4.1ip the initial point of the earlier null geodesic will 
have the coordinates {x,y) — {P,Q), and, at this point, 
£,x — £ ,y — 0. In addition, the isotropy subgroup of G3, 
existing in each case at every point of the manifold, allows 
us to rotate the (x, y) coordinates, with no loss of gen- 
erality, so that the initial value of dy/dr for our chosen 
geodesic is zero, i.e. so that the ray is initially tangent 
to the y = constant subspace. Equation (|4.1ip shows 
that with such initial conditions (and with £,y = at 
the initial point) we have d^y/dr^ = initially, and so 
d'^y/dr^ = = dy/dr all along the geodesic. 

With coordinates chosen in such a way, equations 



For illustration, we first consider a configuration 
that is not realistic, but shows the non-RLP effect in 
a clearly visible way. It is an LT model specified 
by the following functions: ts = and p{to,r) = 
Po [1 + 6 ~ (5exp (— r^/cr^j] , (where to is the current in- 
stant, r is defined as R{to,r), and po is the density at 
the origin and equals 0.3 x (37Jg)/(87rG), where Hq — 
72 km s^^ Mpc"-'^ and G is the gravitational constant). 
This model is the so-called giant void model discussed in 
detail in [23], with the best-fit parameters: S = 4.05 and 
a = 2.96 Gpc. We use this model to study the configu- 
ration presented in Fig. [U where, for the middle curve, 
the angle between the radial direction and the incoming 
geodesic is 7 = 0.227r. We consider 3 light paths. The 
first one corresponds to photons received by the observer 
5 X 10® years ago, the second one corresponds to photons 
received at the current instant, and the third one corre- 
sponds to photons which will be received in 5 x 10® year 
in the future. Figure [1] shows these 3 geodesies projected 
on the space t = now along the flow lines of the matter 



^^ The null geodesies with dx/dr = dy/dr = can properly be 
called radial only in the L— T model, where e = -|-1. What this 
condition means in the other two cases is not clear, so the term 
"radial" is used here only as a brief label. 



o 

Q. 
O 




0.2 0.4 0.6 0.8 1 
[rad] 

FIG. 1: Three nonradial null geodesies in an L-T model, pro- 
jected on the space t — now along the flow lines of the L-T 
dust. Each geodesic runs between the same observer and the 
same emitter, which, at present, lie at the same distance of 
3.5 Gpc from the center. The solid line represents the ray that 
the observer receives at the current instant, 7 = 0.227r, the 
dashed line represents the ray that was received 5 x 10^ years 
ago, 7 = 0.2287r, and the dotted line represents the ray that 
will be received 5 x 10^ years in the future, 7 = 0.2147r. As 
seen, nonradial null geodesies in the L-T model do not have 
the RLP property. The inset shows the density profile (po 
is the density at the origin) evaluated at the current instant 
along these three different paths. 



source in the L-T model. Since in each case the light 
paths are different, the profile of matter density along 
each projected light ray is different. This feature is pre- 
sented in the inset in Fig. [1] Even though the density 
variation along the light path is of small amplitude, the 
effect is clearly visible. The average rate of change of the 
position of the source in the sky, seen by the observer, is 
~ 10~^ arc sec per year. 

Now we will study a more realistic configuration. The 
parameters of the L-T model will be the same as above, 
but the placement of the observer and of the source will 
be different, see Fig. [5J The observer (O) is located at 
Ro [Ro — Rito^r) is the present-day areal distance) and 
observes a galaxy (*), the angle between the direction 
towards the galaxy and towards the origin is 7. We study 
3 configurations: (1) Rq — 3 Gpc, (2) i?o = 1 Gpc, (3) 
Rq — 1 Gpc but with 5 — 10. All 3 cases have d — I Gyr 
(w 306.6 Mpc). For each case (for a given 7) we find a 
null geodesic that joins the observer and the galaxy. We 
then calculate the rate of change 7, which is equivalent 
to the change of the position of the galaxy in the sky. 
A detailed description of the algorithm is presented in 
Appendix [H] The results are presented in Fig. [3l 

As seen, the rate of change of the position of the source 
in the sky depends on the angle 7. The amplitude of 
the change is of the order ~ 10~^ arc sec per year for 
case (2) and ^ 10~^ arc sec per year for cases (1) and 




FIG. 2: A schematic view of the considered configurations. 
The observer (O) is located at Rq and observes a galaxy (*), 
the angle between the direction towards the galaxy and to- 
wards the origin is 7. Because of the non-RLP effect the angle 
at which the galaxy is observed at some other time instant is 
different. 



(3). Given Gaia^"'^ accuracy of position measurement, 
5 — 20 X 10~^ arc sec, we would need to wait at least a 
few years to detect the change of position due to non-RLP 
effects. However, this estimate assumes that we have a 
reference direction that does not change. This will be 
a difficult practical problem, since cosmological observa- 
tions are done under the assumption that our Universe 
is precisely represented in large scales by the Robert- 
son - Walker class of models, in which there is no such 
drift. We would have to identify a direction that does 
not change with time even in an inhomogeneous model 
or measure a relative change of position between various 
objects. 



IX. SUMMARY 

By a method analogous to that of Bondi [ifll we have 
derived the equations to be obeyed by the redshift in a 
general Szekeres P' ^ spacetime, (|5.12p - (|5.14p . The 
null geodesic equations parametrised by r, which must be 
solved together with (|5.12p - (|5.14p . are given by (|i^ - 
(|4.1ip . Although the physically most interesting quantity 
is the longitudinal redshift determined by r, the other 
two components, C and Tp, must be numerically moni- 
tored along the ray because the equations that determine 
(r, C, 'ip) are coupled. 



http://sci.esa.int/science-e/www/area/index.cfm?fareaid=26 
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FIG. 3: The rate of change of position in the sky (7) due 
to the non-RLP effect, expressed as a change of an angle in 
arc sec per year xlO^. The solid line presents case (1) where 
i?o = 3 Gpc, the dashed line presents case (2) where Ro ~ 1 
Gpc, and the dotted line presents case (3) where Ro = 1 Gpc 
and 5 = 10. 



We have shown that, in general, two light rays sent 
from the same source at different times to the same ob- 
server do not proceed through the same succession of 
intermediate matter particles; we refer to this property 
by saying that the Ught paths are not repeatable. In a 
toy model, with the present spatial distance between the 
light source B and the observer being of the order of 1.5 
Gpc, the estimated rate of the drift of B across the sky 
would be ~ 7 X 10^^ arc sec per year. In a more realistic 
configuration, this number is ~ 10~^ arc sec per year. 
The Gaia is expected to have the precision of position 
determination 5 — 20 x 10~^ arc sec. 

We have derived the equations defining repeatable light 
paths ( RLPs) , (lO)) - (|6!^ : they must hold together with 
(1121) - (FTTI) and ((5T^ . We have shown that all null 
geodesies are RLPs only in the Friedmann models. The 
only other cases in which RLPs exist are the following: 

(i) The axially symmetric Szekeres models, in which 
the RLPs are the null geodesies intersecting every space 
of constant time on the axis of symmetry. 

(m) The radial null geodesies in the G3/S2 subcases 
(i.e. in the spacetimes that have 3-dimensional symmetry 
groups). 



Appendix A: Equations of null geodesies in a 
Szekeres spacetime in an aflfine parametrisation 
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Appendix B: Solutions of (|6.5[> . 

Since (|6.5p should hold along certain null geodesies, its 
derivative by r along those geodesies must be zero. This 
derivative, denoted by V/dr, of any quantity x defined 
along the geodesic, x(i(^), ^, a^(''),2/('')), is: 



(Bl) 



Z?Y 


dx dt dx dx dx dx dy 






dr 


dt dr dr dx dr dy dr 



Calculating V/dr of (|6.5p we get: 



For convenience of the readers, the equations of null 
geodesies in a Szekeres spacetime are copied here from 
R-ef. [12| . They are given in an afSne parametrisation. 
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The expression in the last line can be calculated from 
Km - gUl) using dm; it is: 



El 



dh 



E, 



d^x 






dxV 



.(B3) 



Substituting dBS]) and (lO)) - (|i3|) in (|B2]), and taking 



into account the identities S,xy = S,xx ~^ :yy = we get: 

(B4) 



-C/12,r-3— 



da; 
dr 



0. 



This should hold simultaneously with (|6.5p . Since we 
assumed da;/dr ^ ^ dy/dr, (|6.5p and (|B4|) imply: 



E^iE 



12-C/13,r 



£^13-E'12,r — 0. 



(B5) 



When ()44|1 - ()43|) are substituted in f|B5)l . £' factors out, 
and the other factor is: 



^ ir t,^ ly ^ irrx ^ ix ^ -jrry ) 
C ^rr \^ tx ^ ;ry ^^O ,.y C ^j-x } 



^ l^<^ ;rj; ^ irry <-^ iry <^ ^rrx 



) = o. 



(B6) 



This simplifies to a polynomial of second degree in x and 
y, which should vanish identically^^ Using the algebraic 
program Ortocartan [26|, l27| we find that the coefficient 

of {x"^ + y^) is 



r^rr ^ ir ^:r ' 



0. 



(B7) 



One of the solutions of this is P,r — 0; then no limitation 
for Q follows. This case we consider separately below. 
When P,r ^ 0, (|B7| implies 



Q^CoP + Do, 



(B8) 



where Co and Do are arbitrary constants. When this is 
substituted in (jB6[) . the coefficient of y implies: 

£ [SS,r Pirr ^S,r''P,r —SS,rr Pyr ) ^ (1 + Co ) P,r = 0, 

(B9) 
and this guarantees that the whole of (jB6[) is fulfilled. 

The case £ = is seen to be incompatible with P,ry^ 0. 
This means that no RLPs exist in the £ = models with 
P,r 7^ 0. Further calculations apply only to £ = ±1. 

In integrating (|B9I) we can assume 5*,^ ^ because 
S,r = immediately implies P,r = 0, which we have left 
for a separate investigation. Therefore we can introduce 



^^ Should there exist any point Vn in the spacetime at which the 
polynomial would be nonzero, this would mean that the deter- 
minant of the set -f 116. 51 1. | |B4| |1- is nonzero at "Pn, which in turn 
would mean dx/dr = dy/dr = at Vn - contrary to our initial 
assumption. 



S{r) as the new independent variable in (jB9l) . which then 
becomes: 



e {SP,ss -P,s ) - (1 + Co') P,s' = 0. 



(BIO) 



Since the case P = constant was left for later, we assume 
P,s y^ 0, and then (jBlOp is easily integrated with the 
result: 



eS' 



[I + Co') P' = C3P + D-s 



(Bll) 



where C3 and D3 are new arbitrary constants. 

When £ = +1, eqs. (jBSp and (jBlip are equivalent to 
those that were shown in Ref. [12] (sec. 3.3.1) to be 
sufficient conditions for the Szekeres metric to be axially 
symmetric. However, this equivalence is nontrivial, and 
the extension of the proof to £ = 0, —1 is not automatic, 
so we have to elaborate on this subject. 

For this purpose, we note the following properties of 
the general Szekeres metrics (12.71) : 

(1) The metric (|2.7p does not change in form under the 
coordinate transformation: 



{x,y) = {x' + xo,y' + yo), 



(B12) 



where {xo,yo) are arbitrary constants. This changes 

(p,g) to 



iP,Q)^{P-xo,Q-yo). 



(B13) 



(2) The metric (|2.7p does not change in form when 
(x, y) are transformed by a general orthogonal transfor- 
mation: 



ax' + by' —bx' + ay' 

X = , ^ _^ , y - 



Va2 + 62 
which implies the change of (P, Q) to 



Va2 -t- 62 ' 



P 



aP-bQ 

Va2 + 62 ' 



Q 



bP + aQ 
Va^ + 52 ■ 



(B14) 



(B15) 



(3) The metric (|2.7p does not change in form under the 
discrete transformations: 

{x, y) = [y', x'), (x, y) = (-x', y'), 

{x,y) = (x',-y'), (B16) 

which induce, respectively 

{P,QJ^{Q,P), iP,Q)^i-P,Q), 

iP,Q) = iP,-Q). (B17) 

(4) The metric (|2.7p does not change in form when 
(x, y) are transformed by a conformal symmetry of a Eu- 
clidean 2-plane - a 2-dimensional Haantjes transforma- 
tion by the terminology of Ref. [ll| . It has the form: 



x' + Ai (x'^ + y'^) 



T 
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y' + X2(x'^ + y'^) 



T 



(B18) 



dot 



1 + 2\ix' + 2A2?/ + (Ai^ + Aa") (x'^ + /) 



where Ai and A2 are arbitrary constants - the group pa- 
rameters. This group is Abehan, the inverse transforma- 
tion to (IBISP being of the same form, but with parame- 
ters (— Ai,— A2). The characteristic properties of (JB18I) . 
useful in calculations, are: 



da;" 



+ r 



T ' 
dx'^ + dy'^ 

r2 



(B19) 



Under (|BT8)) - ((BT9| . iP,Q,S) change, respectively, to: 



P=-[P-Xi{P^ + Q" 

Q = ^ [g - A2 (P' + Q" 

S = S/U, 

C/ =' I-2A1P-2A2Q 
+ (Ai' + A2')(P' 



-^s')]^ 
-^s')]^ 



Q' + eS^) 



(B20) 



Let £ denote £ with {x,y,P,Q,S) replaced by 
{x' , y', P, Q, S). Then calculation shows that 



£ = £/T, 



(B21) 



and since T does not depend on r, it follows that the Qrr 
component in (|2.7|) is also covariant with (JBISP - (jB19|) . 

Now we will use the properties listed above to interpret 
the consequences of (jBSp and (jBlip for the metric (|2.7p . 

The Do in (|B8| can be set to zero by (|BT21) with 
(a^Oj2/o) = (0,£'o)- The Co in (jB8|) can be set to zero 
by (jB15|) with b = — aCo; the result of these two trans- 
formations is Q — 0- Finally, the C3 in (JB11[) (with 
Co = taken into account) can be set to zero by 
(|B12p with {xo,yo) = (^C3/2,0). Thus we can assume 
Do = Co = C3 = with no loss of generality. 

We carry out a combination of (jB12p with (JBISP : 



X = Xq 
,1 



x' + Ai (x'^ + v'^) 



T 



y = y, (B22) 

and get the following generalisation of (IB20p with A2 = 0: 
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Using (IB8| and ([BTT]) with Do 
becomes 
1 



Cq — C3 



- 0, the above 



P^-[P-xo-Xi {-2xoP + D3 + xo^)] , 



Q = 0, 



(B24) 



dcf 



U = I- 2Ai(P - xo) + Xi {-2xoP + D3+ xo') . 

Now it can be seen that if the constants {xq, Ai), so far 
arbitrary, obey: 



l + 2Ai2:o = 0, 

xo + Ai (L»3 + xo^) = 0, 



(B25) 



then P — Q — 0, and in the {x\y') coordinates the 
Szekeres metric is explicitly axially symmetric. However, 
two things must be noted: 

(1) The set (|B25p has no solutions when D3 < 0. 

(2) With P = Q = 0, eq. ^ is fulfilled identically, 
and (jBlip no longer follows, thus there is no limitation 
onS*. 

Looking at (jBlip with Co = C3 = we see that D3 < Q 
cannot occur when e = +1 or e = 0. The case D3 — 0, 
although possible with e = -|-lor£ = 0, need not be 
considered with these two values of e, for the following 
reasons: With e = +1 this would imply S = 0, which 
is an impossibility in (j2.7p . and with e = = Cq = C3, 
(jBlip implies P = 0. With Q — Q now being considered, 
P ~ e — guarantees that S may be set to 1 by a suit- 
able reparametrisation of the other metric functions [20| . 
Consequently, with P = Q — 0, the e — Szekeres met- 
ric is already plane symmetric even with non-constant S, 
and the Szekeres metrics with 3-dimensional symmetry 
groups are considered in Sec. IVIIl 

So, finally, D3 < must be considered only for e = —1. 
Since these calculations are lengthy and very compli- 
cated, we have moved them to the separate appendix 

We now come back to jBTj to consider the case P,r~ 0. 
By a transformation of x this can be reduced to P = 0. 
Then, the whole of (IB6I) becomes: 



(B26) 
This is equivalent to the subcase Co = of (|B9p under the 
coordinate transformation (x, y) — {y, x) and the associ- 
ated renaming {P,Q) = {Q,P). This case was included 
in the consideration above. 

Thus, apart from the special cases D3 < to be con- 
sidered further on, RLPs with dx/dr ^ 7^ dy/dr may 
possibly exist only when the Szekeres metric is reducible, 
by a coordinate transformation, to one with P = Q = 0. 
In this case, (|6.5I) becomes: 



(B27) 



S,r ( dy dx\ 



U 



But with e = and P — Q — Q now being consid- 
ered, the quasi-plane Szekeres metric is plane symmetric 
even with non-constant S; see the paragraph following 
(jB25p . The Szekeres metrics with 3-dimensional symme- 
try groups are considered in Sec. IVIIl so we need not 
consider e = here. 
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When S,r = = P,r = Q,r, all Szekeres metrics acquire 
a S-dimensional syrametry group and are considered in 
Sec. IVIII Thus, we need not consider 5*,^ = in (jB27[) . 

What remains of (|B27p is xdy/dr — ydx/dr — 0. One 
solution of this is a: = along the null geodesic. The 
other solution is y — Gqx along the geodesic, where Go is 
a constant. However, we are now considering the axially 
symmetric Szekeres solutions in which P = Q = 0. In 
this case, a rotation (|B14|) can be chosen so that y = Gqx 
is transformed to x' = 0. So, apart from the special case 
e = —1,1)3 5: 0, for the other Szekeres solutions the 
following result applies: 

Corollary 2: 

The Szekeres spacetimes in which RLPs exist with 
dx/dr 7^ 7^ dy/dr either are the Friedmann models 
(in which all null geodesies are RLPs) or are inhomo- 
geneous and axially symmetric or have a 3-dimensional 
symmetry group. In the first two cases, coordinates may 
be chosen so that x = along the hypothetic RLP and 
P = Q = m the metric. The second case is considered 
in Appendix If] The third case is considered in Sec. IVIII 
and in Appendix [Gl 



Appendix C: The special metric with Q — 0, e — —1 

and D3 < 0. 

We consider here the special case D3 < that arose in 
solving ()B25p . The Szekeres model in question has 



£ 



2Px + y^ + D3 
25* 



P, 



^12 = j^iy'-^' + Ds) 



E 
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xy. 



(CI) 



The solution of (|6.5p is then either P,r — 0, which belongs 
to the axially symmetric case considered in appendix IE} 
or 



y^-Cy-D3 = 0, 



(C2) 



where C is the arbitrary constant that arises while inte- 
grating (|6.5p . By writing the above as x'^ + {y — C/2)2 — 
D3 + C^/4 we note that the following must hold: 



D3 + CV4 > 0. 



(C3) 



(With this quantity being negative, (jC2p has no solu- 
tions, i.e. there are no RLPs. When it is zero, the only 
solution of (jC2l) is {x, y) = (0, C/2), what is possible only 
in the axially symmetric case of appendix lEl) 

Note that ^^ implies C 7^ 0, since D3 < 0. 

Taking the second derivative of (|C2p by r and substi- 
tuting in it the expressions for x,rr and y,rr from (|4.10p 
- (|4.1ip . we obtain an identity. This means that (|C2p is 
consistent with the geodesic equations (|4.9p - (|4.1ip and 
defines a special class of null geodesies. We will verify in 
the following that this class does not contain any RLPs. 



We note the following auxiliary formulae. In Eqs. (IC4p 
~ (|Clip asterisks mark those equations that hold only 
along the null geodesies obeying (|C2p . those without the 
asterisk are general. After (|Clip all further calculations 
are done only along these geodesies, so we omit the as- 
terisks for better readability. 



f 



(C4) 



-2Px + Cy -I- 2D3 
2S 

X — r y Xr^r ^^r p 



s ' '" s 



s s 

(C5) 
..--%, (C6) 
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P, 



W ^12 = TT^^ (2y' - Cy) 



25*2 



(*) y,r^- 



2x 



2y-C 



X,r , 



^ ' ''"■ ^^'^ ~ {2y-CY ''"■ • 

(*) 2x,r'^Ei2+2x,rV,rEi'i 

^^y(2y~C){x,r^ + y,r^) 



From (|5.15p we have 



^•>r \ y ir ^ir 



•^1 



£^ 



$2 (e-fc)$2- 



(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 



In order to write the equations in a compact way it 
is convenient to use the symbol ^ introduced in (|6.4p . 
Recall: this is a factor of shear, and when it vanishes, the 
Szekeres model reduces to Friedmann. Thus, in searching 
for RLPs in nontrivial models we will assume ^ 7^ 0. 
Using (jClip in (|5.16p adapted to RLPs we obtain 



"^ir ^,r — 



We also have: 



£-k $ {e-k)^ 

e - k $ 



(C12) 



$01 - $,t *l/$ = $,tr -$,t $,r /$ = *, (C13) 
$,«r-$,U$,r/$ = *,t+^*, (CM) 

$,t„-$,t$,„./$-*,,+^^. (C15) 

Assuming vli 7^ wc now multiply (16. 2p by $$i^i,r /[(e — 
A:)*], use (|C7p - (jClSp . cancel r that factors out, and 
write the result in the form: 

X,r {t,/ + C2t,r^ + 03*,^ +C4) = Bit J + B^t.r , (C16) 

where 

dcf 2$($i*,t-*2) 



C2 



(e-fc)* 



(C17) 
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C3 

C4 

Bi 
B3 



dcf -3$!^ + ^^i^„. /I' - $$,r 

e — k 



e — k 



dcf 2$$!^^ 

(e - fc)2 ' 
dcf P,ry{2y~C) 



dcf 



(e-fc)S'2 
'2{e~k) 



Bi. 



(C18) 
(C19) 
(C20) 
(C21) 



Then, using jCll, jCTl) and (jClip we can rewrite (IT^ 
in the form: 



where 
C5 

C7 



C5i,r^ + cetjr^ + c^t^r +C8 + Aa;,^ t,r , (C22) 



dcf 



dcf 



dcf 



dcf 



dcf 



$$1 ' 



$1 



$, 






$1 

(47^3 + C2)P„y$ 



2(e-fc) 



(2y-C)S'2£2$^ • 
Combining ((C9l) and (ICllI) we get: 



^ir 



(2;/-C)^g^ 

(4£l3+C2)$2 



^ir 



e ~ k 



(C23) 
(C24) 

(C25) 
(C26) 

(C27) 



(C28) 



Equations (|C16|) and (|C28p determine di/dr along the 
hypothetic RLP. Formally, a solution for di/dr of these 
equations always exists, but it must be consistent with 
the geodesic equations, and this is what we will inves- 
tigate next. Namely, every solution of these equations 
must be preserved along the null geodesies. To see 
whether it is, we first transform this set into a single 
polynomial equation for dt/dr. 

We square (jC16p and use (|C28I) in the result. We thus 
obtain an 8-th degree polynomial in i,r, whose coefficient 
atj L,^ IS 



2c2 



Ckl = 



(2y~C)2g 
(4D3-f C2)$2- 



(C29) 



It is seen that it cannot vanish except when y — C/2, but 
this defines a "radial" geodesic that exists only in the ax- 
ially symmetric case [l2| . Thus we divide the 8-th degree 
polynomial by ai and obtain the following equation 



Q '7 fi ^ A "^ 

t,r + '2c2t,r + 0-3^, r + a^t^r -f a^t,r -\- O-et^r 



+a7t,r + ast,r -I- ag = 0, 



(C30) 



where: 

dcf , 2 '^1' (4i?3 + C2) $2p„,2y2 

03 = 2C3 -I- C2 



e-k (e - /c)2S'4f 2 



dcf „ , „ 0^1 



£ — fc' 



(C31) 
(C32) 



dcf r, , 2 ^r. , 2\ *^1 

as = 2C2C4 + C3 ~ [2c3 + C2 * 



e — k 



-3^lL<i5i±SS^r!»!. (C33, 



(£ - fc)3 



S'4f2 



ag = 2C3C4 - (2c4 + 2C2C3) -, 



e — k' 



clct 2 /^ 

ay = C4 - (2C2C4 -f C3 



$1^ 



e — k 
91>i" (4J3 + C^)$^P,.V 

4(£ - fc)4 S'4£2 

2 



dcf „ *1 

= -2C3C4 



dcf 2 *^1 

ag = -C4 



e — k 
2 



e — k 



(C34) 

(C35) 
(C36) 

(C37) 



Now we differentiate (IC30I) along the null geodesic by the 
rule (jBll) . and use (|C30|) to eliminate <,r^°, i,r^ and t,r^ 
from the result. In this way we obtain: 

bit, J. -\- 02t,r + "S^jr ~l~ b^tjr + O^tjr "I" "6^jr 

-f a;,^ (;9lt,/ + /32i,r^ + /33i,r^ + /?4t,r^ + /35i,r^ 

-|-/36t,r' + P7t,r + /Ss) = 0, (C38) 

where 

dcf 

61 = 8C8 -t- 2c2,r + 03,4 — 2a3C6 — 3a4C5 — 2C2C7 

-4c2C2,t -I- 6a3C2C5 -f 4C2^C6 - 8C2^C5, (C39) 
def 

62 = 03, r + ai,t + 14c2C8 — 2a3C7 — 3a4C6 — 4a5C5 

— 2a3C2,t -I- 2a4C2C5 -I- 2a3C2C6 -I- 203205 



-4a3C2 C5, 



(C40) 



dcf 



63 = a4,r + as,* + 6a3C8 - 3a4C7 - 4a5C6 - haec^ 
—2a4C2.t + 2a5C2C5 + 2a4C2Ce + 2030405 



-4a4C2 C5, 



(C41) 



dcf 



64 = 05, r -t- 06,4 + 5a4C8 — 4a5C7 — bagCQ — 60705 

— 205C2,t + 2O6C2C5 -I- 2O5C2C6 -I- 2O3O5C5 



-405C2 C5, 



(C42) 



dcf 



&5 = a6,r + a?,* + 405C8 — 5o6C7 — 607C6 — 7o8C5 
— 206C2,t + 207C2C5 -I- 2O6C2C6 + 2O3O6C5 



-406C2 C5, 



(C43) 



dcf 



&6 = <^7,r + ^8,* + 3ogC8 — 607C7 — 7o8C6 — 8ogC5 

—2arC2.t + 208C2C5 + 207C2C6 + 20307C5 
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-4a7C2 C5, 



(C44) 



dcf 



7 LlUl 

O7 = 



os.r + cig.t + 2a7C8 - 7asCj - SagCe ~ 2a8C2,t 

+ 2agC2C5 + 2a8C2C6 + 2030305 - 4o8C2^C5, 



dof 



&8 = a9,r + asCS - SogCT - 209C2,t + 2O9C2C6 
+ 2030905 - 409C2 C5, 

A == 2 ( C2,a; - _ „ C2,y j ~ 2C2A, 

o dof 2x 

P2 = a3,x - _ „ a3^.y - 2o3^, 

o dof 2a; 

P3 = a4,a; - T^ J^ a^.y - OO4A, 

zy — O 

« dof 2a; 

P4 = 05^2: - _ 05^,y - 405^4, 

« dof 2x 

o def 2x 

P6 = a7,x - — — -p^ aj^y - bajA, 



« dof 

P7 = 



2y~C 
2x 



as^y - 7asA, 



o dof 2a; 

P8 = ag^x - „ _ „ a9^.y - 809A. 



(C45) 

(C46) 

(C47) 

(C48) 
(C49) 
(C50) 
(C51) 
(C52) 
(C53) 
(C54) 



We provisionally assume that the coefBcient of a;,^ in 
(|C16p is nonzero. (We will later come back to this point 
and investigate what happens when it is zero.) Then 
we determine x,r from ()C16p and substitute the result 
in (|C38p . After multiplying out to get a polynomial in 
i,r, we again use (|C30I) to eliminate i,r^° and t,r'^ (but 
not t,r^)- Then we assume that the coefficient of t,^^, 
denoted di, is nonzero (we will check the case di = 
later), and divide the equation by di. In this way we 
obtain: 

+ S7t,r^ + S8t,r +S9 = 0, (C55) 

where Si = di/di, i = 2, . . . , 9, and 

dof 

di = 63 - &2C2 + &1C3 + /3ii?3 + l^sBi - ash 

-asPiBi - 2/32B1C2 + 2biC2^ + 4^lBlC2^ 

(C56) 

dof 

^2 = bi + hiCi + 63C2 + 62C3 + ^2-83 + jS^Bi - 0461 
-0362 - aiPiBi - a-il32Bi + 0361C2 
+2a3/3iBiC2, (C57) 

dof 

^3 = ^5 + hci + 64C2 + &3C3 + l3zB'i + jS^Bi - 0561 
-0462 - a^sPiBi - 04/32-81 + 046102 
+2a4^iBiC2, (CSS) 

dof 

di = 65 + &3C4 + 65C2 + &4C3 + PaBs + PeBi - aebi 
-0562 - aePiBi - 05/32-81 + 0561C2 
+2o5/3iBiC2, (C59) 



dof 

^5 = bj + biCi + b6C2 + &5C3 + ^^5-83 + ^^7-81 - 0761 
-0662 - a7;Si-8i - 06/32i?i + O661C2 
+206A-81C2, (C60) 

dof 

de = bs + 65C4 + &7C2 + beC3 + jSqB^ + l3%Bi - 0361 
— 0762 — a-sPiBi — 07/32 i?i + 07&1C2 
+2o7/3iBiC2, (C61) 

dcf 

dj = bgCi + &8C2 + ^723 + ;57-83 - 0961 - 0362 

-ag/3iBi ~ asl32Bi + 036102 + 2o8/3ii?iC2, 

(C62) 

dcf 
^8 = &7C4 + &8C3 + /38-83 - O962 - 09/32-81 

+0961C2 + 2o9/3ii?iC2, (C63) 

(C64) 



dg = 63 C4 



Every solution of (JC30I) is a candidate RLP, and every 
RLP must obey (jC30p . Equation (jC55[) is the condition 
that (|C30p is preserved along null geodesies. Thus, every 
solution of (|C30|) must also be a solution of (jCSSp . Since 
(|C30p and (jC55|) are of the same degree in t,r, it follows 
that both must have the same set of zeros. Consequently, 
their coefficients must be the same. After we make sure 
that they are the same, we may next investigate which 
zeros define RLPs. Thus, the following equations are 
necessary conditions for the existence of RLPs: 



2C2 - (52, 

-4=4> 2c2(ii 



(Ji'i f^t 1 ' "J; ' 

d2 — 0, diUi - 



..9, 
d, = 0. (C65) 



By far the simplest condition, as seen from (jC64p . is the 
one with i = 9. Even so, further calculations are so 
complicated and involve intermediate equations so large 
that they could be done only using the computer algebra 
system Ortocartan [26l . \2l\ . and we only describe how 
they were done. 

First we observe that the functions ($1,2/, f) are lin- 
early independent, and in the resulting final equation can 
be used as independent variables. Although the proof is 
a simple exercise, it requires careful inspection of special 
cases that we had earlier excluded for separate investiga- 
tion, so we give it in the separate Appendix |D] 

The condition (|C65P corresponding to z = 9 is 



diog — ^9 = 0. 



(C66) 



In this, one has to do the whole cascade of substitutions, 
listed in ((CTT)) - ((C64)) . In the result, we use (IC5|) - ([Ce]) . 
However, we use the last of (jCSp only to eliminate £,„■ 
For £,r we substitute from (|4.ip . i.e. 



£•,, = £($,,-$1)/$, 



(C67) 



in order to express £,r through $1. 

We then use (JCISP to express $,trr through VP, and 
(16. 4p to express $,tr through 'if. In the result we use 
(IC2p to eliminate x^ and (jC4l) to express x through £. 
From the resulting equation we can factor out $1 , and 
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we must multiply it by £^ to get rid of negative powers 
of £. The final equation thus obtained has on its l.h.s. 
a polynomial of 4th degree in £, of 4th degree in y and 
of 6th degree in <I>i (recall, we determined that ($i,y,£) 
are independent variables). ^^ So, if this polynomial is to 
be zero, the coefficients of all powers of the independent 
variables must vanish separately. 

Now we take this large polynomial as data for a second 
program, in which (^i,?/,,?) are treated as independent 
variables, no longer as functions. In it, we determine the 
coefficient of y"^ and substitute £ — Q to find the term 
independent of £. The resulting equation is: 

16 (4L>3 + C^f $^**,t P,// [S^{e - kf] = 0. (C68) 

In this equation we can discard several alternatives: "I> = 
obviously, AD^+C'^ = because of (|C3)) . * = because 
it defines the Friedmann limit and P,r — because, in 
view of (JBIII) and the paragraph above (jB22p , it leads to 
the G3/S2 symmetric cases considered in Appendix [G] 
The only case to consider is thus ^,f = 0. 

In order to investigate it, we substitute ^,t = in the 
large main polynomial, and in the resulting somewhat 
smaller polynomial we take the coefficient of j/^. The 
equation that results is: 

-16 {4D3 + C^) (4P2 + c^^ £$5^3 

xPj/ [PS\e - kf] = 0. (C69) 

The factors £ and $ obviously cannot vanish, and why 
zero values of (4D3 + C^), ^I' and P,r are discarded was 
explained above. Thus, the only case left is 4P^+C^ = 0, 
which means P = = C. But this is just a special case of 
P,r ~ discarded above. Thus, (|C68p does not include 
any case that would define any new RLP, apart from 
those considered elsewhere. 

We go back to (|C55P to consider the case di = 0. The 
calculation is almost the same as we did for (|C66I) . with 
only minor differences: this time the expression is some- 
what simpler, and $1^ does not factor out. We employ 
the algebraic program to calculate £^di, with the same 
cascade of substitutions as before, take the coefficient of 
y'^ at £ = 0, and obtain an equation almost identical to 
(IC681) : 

4 {4D3 + C'^f $4*,t F,// [5«(e - fc)%] = 0. (C70) 

As explained above, only *,t could possibly be zero, so 
we substitute ^,f = in the main large polynomial, and 
in the resulting expression take the coefficient of y"^ . The 
result is almost the same as (IC69I) 

-4 [AD^i + C^) (4P2 + (^2) £$3^ 



and again does not include any case that would define a 

new RLP. 

Finally, we go back to (|C54|) . where we assumed that 
the coefficient of x,^ in (|C16p was nonzero, and investi- 
gate what happens when it is zero. Then 



Bit^r + B^t 



S^ir 



0, 



(C72) 



and one of the solutions of this is i,r = 0. This we im- 
mediately discard because it defines a spacelike curve, 
while our RLPs must be null geodesies. In consequence 
of (|C2T|l . another solution of ([C72|) is Bi = 0. But this 
implies P,^ =0or2/ = 0orj/ = C/2. The first case 
leads to the G3/S2 solutions considered in Appendix [Gl 
the other two to the axially symmetric solutions of Ap- 
pendix |F1 So the only possibility left to fulfil (|C72p is 



t,r — ±1 



$1 



(C73) 



2{e-k) 
Putting this into the coefficient of x,r in (|C16|) we get: 



<- 


■ 


3/2 r 
$1 < 

23/2$^ 




-k)\ 




3 


v/3(£-fc) 



Tl 



± 



23/2$^^j 

v/3(£-fc)* 



*l^ 



$ 



4* (*¥-*•'■• 



= 0. 



$, 



$1 



(C74) 



In this expression, we do the same series of substitu- 
tions that we did in the large polynomial that resulted 
from (IC66P : we express £,„ through £,r using (|C5p . £,r 
through $,r and $1 using (jC67p . then x through £ and y 
using (jC4p . and multiply the whole expression by £. The 
result is: 



Tl 



23/2$^,t 

v/3(e-fc)* 

23/2$\]/ 



£^i' 



2 / *,r 
± _ $^1 

3^3{s-k) 3 V * 

2 2 o r/p„ 



\ C j_ '*^ '^ 



^S£_ 
"~P" 

= 0. 



Cy 
2P 



£$1 



P 



(C75) 



We then use the fact that ($1, j/,f ) are linearly indepen- 
dent and require that each coefficient of an independent 
function is zero. There is only one term containing y, 
with the coefficient C{P,r /S),r- But C cannot be zero, as 
explained below (IC3p . Thus {P,r /S),r— is the unique 
implication of this (it will be seen from the following that 
we need not consider whether this condition is consistent 
with the other equations that P and S must obey) . This 
means: 



The whole equation would take 1830 print hncs on paper. 



P,r = aoS, 



(C76) 
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where ag is an arbitrary constant. Taking this into ac- 
count, and taking the term independent of $1 in (|C75p 
we get: 



*,. 



-*i^ 



S,r $,r 

5* 



= 0. 



(C77) 



Using (JC76I) . this is easily integrated with the result: 
Xi{t)P , X2{t) 



$ 



aoS 



S 



(C78) 



where {xi{t),X2{t)) are arbitrary functions of t. Both 
appear as integration "constants" of (|C77p . Using such 
$ in the definition of ^, ()6.4|) . we get: 



*$ = 



7(i) 



, s dcf 
lit) = X2Xl.t 



XlX2.t- 



(C79) 



But with ^ = 7(i)/(5$) the last three terms in the coef- 
ficient of $1 in (jC75P sum up to zero, and what remains 
of that coefficient is the equation $^ = 0. The only so- 
lution of this can be ^ = 0, but we know it leads to the 
Friedmann model. 

Consequently, the coefficient of x,r in (|C16p is always 
nonzero. 

Since (|C68|) and (|C70P were, in their respective 
branches of the calculation, among the necessary condi- 
tions for the existence of RLPs, we conclude that the spe- 
cial quasi-hyperbolic Szekeres solution defined by (|C1|) 
does not contain any RLPs except (possibly) when it be- 
comes axially symmetric or G3/S2, but these cases are 
considered in Appendices |F] and [G] 



Appendix D: Proof that {^i,y,£) are linearly 
independent 



We take the equation 



a$i + (3y + •yS ^ 



(Dl) 



with constant coefficients a, /3, 7 and prove that it implies 
a = ^ = 7 = 0. 

We substitute for $1 from (14.11) . then multiply the 
equation by 2S£ and use (|C2p to eliminate x^. We thus 
obtain a polynomial of degree 1 in a; and degree 2 in y, 
which we denote by P. We take the second derivative of 
V hy xy. The result is: 



2P{I3S - 7C) = 0. 



(D2) 



We discard the solution P — because this implies con- 
stant S (see (|Blip and the remarks above ()B22P ). and 
then the metric acquires a G3/S2 symmetry group - these 
are discussed in Appendix [Gl We also discard the case 
(3^0 because then S, and consequently P, are constant, 
again leading to the G3/S2 case. So finally, the implica- 
tion of ini]) is 



^ = = 7C. 



(D3) 



With this, we go back to P and take its second deriva- 
tive by y. The result is —A-fP'^/S = 0, and the solution 
of this is 7 = 0. We again go back to P with /3 = 7 = 0, 
and take its first derivative by y. The result is: 



aC{<^,r+'^S,r/S)^0. 



(D4) 



When the expression in parentheses vanishes, $ becomes 
a product of the form R{t)/S{r), where R{t) is an ar- 
bitrary function. Such form of $ defines the Friedmann 
limit, in which we know that all null geodesies are RLPs, 
so we discard this case. Thus, we follow the case aC = 0. 
Putting this, together with (3 = j — 0, inP and taking 
the derivative of the result by x we obtain 



2a {-P^,r +Pr $ - ^PS,r / S) = 0. 



(D5) 



If the expression in parentheses should vanish, then the 
solution is $ = R{t)P{r)/S{r), which again leads back to 
the Friedmann model. Thus, finally, ()D5[) implies a = 0, 
which completes the proof. 



Appendix E: The RLPs with P 
the geodesic. 



and X = along 



We can leave aside the case when Q,r = because then 
the metric is axially symmetric from the beginning. 

It is useful to turn this case back to that of Appendix 
IH by the transformation (IBM)) - (|BT5)) with a ^ j^ 
b. After the transformation we have x' = —hy' /a, i.e. 
dx'/Ar 7^ if dy'/dr ^ O^and P = -6Q/a, i.e. P,r 7^ if 
(9,.r 7^ 0. Thus the new P and Q obey (HH with -Do = 
from the beginning, and the RLP condition reduces to 
(iBTTj) alone, with (P, Q) replaced by (P, Q). The rest of 
the reasoning of Appendix [Bj then applies, unchanged, to 
(P, (5), with the same result: 

Corollary 3: 

RLPs with P — Q and x = along the geodesic may 
exist only in the special case when the coordinates may 
be transformed so that Q = as well, i.e. the metric 
is axially symmetric, or has a 3-dimensional symmetry 
group. 



Appendix F: The axially symmetric case P = Q = G: 
only the axial geodesies x,r= y,r= are RLPs 

We know from Ref. i]^ that in the quasi-spherical 
case e = -f 1 null geodesies on which x and y are constant 
exist only when the Szekeres model is axially symmetric. 
Then coordinates may be chosen so that P = Q = 0, and 
the constant- (x,2/) null geodesies have x — y — 0, i.e. 
intersect each t = constant space on the symmetry axis. 

In this appendix we show that the statement above ap- 
plies also with £ = and e = —1, that the constant-(a;, y) 
null geodesies are RLPs, and that other RLPs may exist 
only when the Szekeres spacetime has more symmetries. 
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1. Constant-(a;, J/) null geodesies exist only in the 
axially symmetric case 

This thesis was proven in Ref. [l2] for e = +1. The 
assumption made there in the proof, that f > 0, does 
not hold for e = and e = —1, so we first verify what 
happens when £ = 0. 

It is seen from (I4.10[) and (|4.11l) that constant (x, y) 
imply Ei2 — £'13 = along the geodesic. With £ — Q 
(|4.4p and (|4.5|) then imply that either £,r = at aU r, 
which means a G3/S2 symmetry (discussed in Appendix 
[Gt . or £,x— £iy— along the geodesic, which means 
P and Q being constant, i.e. axial symmetry. Thus, 
E — Q along a constant-(a::,2/) null geodesic implies axial 
symmetry anyway. ^'^ 

The equations of Sec. 3.3.1 in Ref. [i2| that are im- 
posed on {P,Q,S) by the condition of constant {x,y) 
along the geodesic become subcases of our (jB8[) and 
(|Blip for a general e. As shown in our Appendix |BJ 
they imply axial symmetry for any e. This is true even 
for the special solution discussed in Appendix [Cl as we 
now show. 

When x,r — y,r— along a nuU geodesic, as stated 
above, (j4.10p and (I4.11|) imply £'12 = £13 = along this 
geodesic. Then (jCip implies that either (i) P,r = 0, or 
(ii) X — and y'^+D^ = 0, or (Hi) y — and x^ — D^ — 0. 
Case (i) is axially symmetric. Case (ii) implies £ = 
along the geodesic, and this was discussed above. Case 
{Hi) implies D3 > 0. However, the solution of Appendix 
[Cl has D3 < by definition. So the only subcase to 



consider here is D3 — : 
But then we have again £ 



^ X — along this geodesic. 
0, which completes the proof. 



2. Constant-(a;, y) null geodesies are RLPs 

As stated above, along null geodesies of constant {x, y) 
we have £12 = £13 = 0. Then (H^) and dOJ are fulfilled 
identically, which means that such geodesies are RLPs. 



3. Other RLPs may exist in the axially symmetrie 
ease only Avith higher symmetries 

We will now show that, in the axially symmetric case, 
(|6.2p - (|6.3|) may have other solutions than constant (x, y) 
only when the spacetime has more symmetries than just 
the axial. 

The whole reasoning and calculation is closely analo- 
gous to the one presented in Appendix [Q for the special 
Szekeres solution. We proved in Appendix |B] that in the 



'^^ Moreover, as shown in Ref. |2C| . the location £ = is infinitely 
far from any point within the spacetime, i.e. does not in fact 
belong to the spacetime. 



axially symmetric case coordinates may be chosen so that 
P — Q — Q and the candidate RLP has a; = 0. Then: 

_ x^+y^ 1 
^ - 25 + 2 ^^' 
£12 = exS,r /S, £13 = eyS,r /S. (Fl) 

Note that with e = 0, this axially symmetric Szekeres 
solution is in fact plane symmetric. Thus, it will be con- 
sidered together with other G3/S2 symmetric solutions 
in Appendix \G\ From here on in the present appendix 
we assume £ 7^ 0, i.e. e = ±1. 

With dni obeyed, (jg ^ an d (p3t are fulfiUed identi- 
cally along X — 0. From (J5.15I) we obtain: 



y,/ ^ hi:___^l__ 

£2 - $2 (e_A:)$2' 



(F2) 



and (jC12[) applies unchanged. We then multiply 
by $$i^t,^ /[(£ - k)^] and use (|F2|), (ICTl - (lUfsl) and 
a; = in the result. As before, r factors out and is 
cancelled, and we obtain an equation almost identical to 
(JC16I) . with the same definitions of (-63,02,03,04), but 
with y,r in place of x,r, and with the definition of Bi 
changed to: 



Bi 



dct 2eyS,r 



{e~k)S' 



(F3) 



We proceed in strict analogy to Appendix |B] From 
(|4.9p using (jF2p we again obtain (jC22p . with y,r in place 
of x,r, and with the same definitions of (05, . . . , cg), but 
with the definition of A changed to: 



^ dcf 2ey<PS,r 



S£^<Pi 



(F4) 



Then we square the current analogue of (jC16p and use 
(|F2p to eliminate j/,r^ from the result. We obtain an 8-th 
degree polynomial in t,r, but this time the coefficient of 

t,/ is 



2 //f,2 



ai = £7$ 



(F5) 



and is sure to be nonzero. Dividing the polynomial by ai 
we obtain (jC38p again, but with the definitions of some 
of the coefficients changed as follows: 



as = 2c3 + C2 



e~k (e-/c)2S'2£2' 



05 = 2C2C4 + C3 — (2C3 



C2 



$1^ 

e — k 



(£ - fc)352f 2 ' 



(F6) 



(F7) 



r-2ff,2ff, 4q 2^,2 



0.7 



dcf 2 to ^ 2^ *1 18£2$^$i^5,/2/2 

= C4 - (2C2C4 -t- C3 j 



e — k 



(£ - fc)452f2 ' 

(F8) 



the remaining ones are the same as given by (jC32p , (|C34p 
and (IC361) - (fC37)) . 
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Now we differentiate the current analogue of (jC30|) by 
r along the null geodesic by the rule (jBip . This time, 
however, a; = along our candidate RLP, so no coeffi- 
cient depends on x. We then use our analogue of (jC30p 
to eliminate i,r^°, t,r^ and t,r^ from the result. The 
equation that emerges is an analogue of (|C38|) with y,r 
in place of x,r, with the same definitions of (&i, . . . , bs), 
and with the definitions of (/3i, . . . , /Sg) changed to 



dcf 



/3i "=^' 2c2.-2c2A 



P2 = as, a - 


- 203^, 


fl dcf 

P3 = 04, a - 


- 804^, 


„ dcf 

Pi = ab,y - 


- 405^, 



dcf 



dcf 



fle,' 



SflgA, 



Pa ~ 0.1, y ~ 607^1, 



dcf 



/?? = fls.a - ^asA, 



/38 



dcf 



fl9,y — oflgA, 



(F9) 
(FIO) 
(Fll) 
(F12) 
(F13) 
(F14) 
(F15) 
(F16) 



where for A the definition (jF4p must be used. 

Here we can assume that the coefficient oi y^r in the 
present analogue of (ICISP is nonzero - the explanation 
given in the paragraphs containing (jC72p - (jC75p still 
applies, except that the B\ given by (jF3l) cannot vanish 
for somewhat different reasons. ^^ Then we determine y,^ 
from that equation and substitute the result in the cur- 
rent analogue of (|C38I) . After multiplying out to get a 
polynomial in i,^, we again use (jC30p to eliminate t^r^^ 
and t,r^ (but not i,r*)- Then we assume that the coef- 
ficient of t,r^, denoted di, is nonzero (we will check the 
case di = later) , and divide the equation by d\ . In this 
way we obtain an exact copy of (jCSSp . with the same def- 
initions (JC56I) - (jC64p of the coefficients; but it is to be 
remembered that some of the symbols in these formulae 
(namely i?i, -B2, 03, 05, 07 and all of (/3i, . . . ,/38)) now 
have different definitions from those in Appendix [CJ 

Consequently, eqs. (|C65I) must still hold, and we again 
choose (|C66P to investigate, by exactly the same method 
as before. By the method of Appendix |D] we show that 
($1, y, £) are still linearly independent in the present case 
(i.e. with £ given by (jFll) . and along the a; = geodesic). 
The explanation given under (IC67I) still applies, with the 
modification that now x is nowhere present, so does not 
have to be eliminated. In place of (jC68p we now obtain: 



256e^$6**,t 5,// [5^(e - kf\ = 0. 



(F17) 



We recall that we excluded the case e = (since it is 
treated in Appendix[G|) , and \1/ = because it reduces the 



^^ The cases £ = and S,r = define metrics of higher symmetry, 
treated in Appendix[G] while y = (together with x = assumed 
throughout this appendix) defines an axial geodesic, which we 
already know is an RLP. 



Szekeres model to Friedmann. We can also exclude 5,^ = 
because then the metric acquires a G^/ S2 symmetry 
and is also treated in Appendix [G] So, as before, the 
only case left to investigate is ^,t — 0. 

We substitute '^,t— in the large main polynomial, 
and in the resulting smaller polynomial we take the co- 
efhcient of y''. The equation that results is: 



2bQe^£<^^-^^S,r^/ [S^{e - kf] = 0. 



(F18) 



The only factors that could vanish here are ^ and 5",^, 
but, as explained above, their vanishing leads to simpler 
cases of higher symmetry. Thus, (jF18p does not include 
any case that would define any new RLP, apart from 
those considered elsewhere. 

We go back to the paragraph after (|F16p to consider 
the case di = 0. The explanation given above (|C70p still 
applies, but this time, in the expression £'*di calculated 
by the algebraic program, we take the coefhcient of y'^ at 
£ — Q, and obtain: 



64e4$%,j 5,,V [S^{e - kf^ = 0. 



(F19) 



As explained above, only ^,t could possibly be zero, so 
we substitute ^,t = in the main large polynomial, and 
in the resulting expression take the coefficient of y*. The 
result is: 



64£3£5$3^5',// [5-4 



kf 



0, 



(F20) 



and again does not include any case that would define a 
new RLP. 

So, the final conclusion is: 

Corollary 4: 

In the axially symmetric Szekeres solutions, apart from 
cases of higher symmetry, the only RLPs are the axial 
null geodesies that intersect each 3-space of constant t 
on the symmetry axis. 



Appendix G: There are no non-radial RLPs in any 
G3/S2 model. 

We will investigate the equation x = (see (|7.ip ') 
and will show that it has no solutions defining nonradial 
RLPs, unless the model reduces to Friedmann. 

Several equations in this Appendix follow from the cor- 
responding ones in the Appendix [C] as the special case 
£,r = 0; they are similar but not identical. 

We will use all equations adapted to the special case 
discussed in Sec. IVIIl i.e. ( = ip — £^— r] = £,r = 
dy/dr = 0. From (|5.15p we find 



£2 



t 2 

"$2" 



<^r 



(e-A:)$2- 



Then, using (jGip in (|5.16p we obtain 



'^;r t^;r " 



$ 



t. 



(Gl) 



(G2) 
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where ^ is defined by (|6.4p - as seen from (I2.9p this is 
a coefficient of shear, whose vanishing defines the Fried- 
mann limit. 

We now substitute (IGll) and (|G2p in x = 0, where x is 
given by (|7.1I) . We multiply the result by $$,r^i,r /[(e - 
fc)^], and cancel r that factors out. The result is: 



Wi =^ t,r^ + C2tJ + C3t,r +C4 = 0, (G3) 



where 



def 2^($„.^,,-^^) 



(G4) 



dcf $($,r*,r-$,rr*)-2$,r* ,^^^ 

C3 = TTTT^ (G5) 



(e-/c)* 



C4 



dcf 



2$$,/* 

(e-A:)2 ■ 



(G6) 



Adapting (|4.9p to the G3/S2 case, eliminating a;,^ with 
use of (IGll) and using (|G3p to eliminate i,^'^ we obtain: 



irr = CQt,r + Cjt^r +C8, 



where 



dcf 



2*,t $, 



C6 

C? = 



dcf *,r $,r , fc,; 



C8 



* $ 2(e-A:)' 

drf ^,r-^ 

£ — k 



(G7) 

(G8) 

(G9) 

(GIO) 



Now we differentiate (|G3p along a null geodesic (since 
the equation must hold all along it), by the rule given 
in (IBip . and use (|G7I) to eliminate t,„- The resulting 
equation is of 4-th degree in i,^. We eliminate the 4th 
power of t,j. by using (IG3I) . In the end, we obtain an 
equation of degree 3 in i,^, which we write symbolically 
as follows: 

dit,r^ +d2t,r^ +d3t,r+di^Q. (Gil) 

The expressions for the coefficients in (jGlip are: 

di = C2.t + 3c7 - C2C6, (G12) 

^2 = C2.r + Cs^t + 3C8 + 2c2C7 - 2c3C6, (G13) 
^3 = C3,r + Ci,t + 2C2C8 + C3C7 - 3C4C6, (G14) 
^4 = C4,r + C3C8. (G15) 

For the beginning, let us assume that di ^ 0. For 
further considerations it will be more convenient to write 
((GTT|) as follows: 

W2 =^ i,r^ + 52t,r^ + (53t„. +,^4 - 0, (G16) 

where 5i — di/di, i — 2,3,4. 

Equation (|G3p is equivalent to x = in (|7.ip . Thus 
(IG3P . just like (|7.ip . defines the collection of RLPs to- 
gether with the conditions of their existence. Every solu- 
tion of (|G3p and (|7.ip is a candidate RLP, and every RLP 



must obey ((G3|) and dLB. Then, (|G16p is the condition 
that the solutions of (|G3p are preserved along the null 
geodesies, thus every solution of (jG3p must be a solution 
of (|G16P and vice versa. But if (|G3]) and (|G16p have 
the same set of solutions, then they must be identical, 
i.e. their respective coefficients must be equal. Thus, the 
necessary conditions for the existence of RLPs are: 



= S,, 



Cidi — rfj = 0, 1 = 2,3, 4. 



(G17) 



As with the previous calculations, we employed the al- 
gebraic program Ortocartan [2y, [23] . We consider (jG17p 
with i = 4, the simplest one. In it, we substitute (jG4p - 
(JGISP and multiply the result by (^/$,,. ) to get a poly- 
nomial in $, ^ and their derivatives (<I>,r factors out in 
the original expression). The resulting expression is sim- 
ple enough to be shown here: 



W.= 



(e - kf 
(e-fc)3 

(£-fc)2 



0. (G18) 



We assume A = and take the fc > model for the 
beginning. The solution of (|2.3p can then be written as 



$(<,r) 
rj — sin 77 



T 

fc3/2 

Af 



(1- cos 77), 



(G19) 



where 77 is a parameter (dependent on t and r), and tsir) 
is an arbitrary function, the bang time. We introduce the 
abbreviations: 



Im 



dcf 



OrZ ^j 



M ' 



D 



def rC Cfj^ 



A/ 



2fc M ' M 

The derivatives of $ and ^I' can then be written as 



(G20) 



*"^= ( "T ) ,r(l-cosr7) 



M/jvf sin 77(77 — sin 77) 



cosry 



MD 



M 



SHUT] 



$, 



fc3/2 2. — cos 1] ' 

vfc sin 77 



1 — cos 7/ 



(G21) 
(G22) 






-I- ( — ) ,^sm?7 



(1 — cost;) 

77 — sin 77 



21. 



D 



M 



M- 



+ 



+ 



1- COS 77 Vk{l-cosr]) 
M{lM),r sin 77(77 - sin 77) 

k 1 — COS T] 

MIm^ (2 sin 77 — sin 77 cos rj — r]){r] — sin 77) 



k 



(1 — COS 77)2 
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MImDm 3 sin -q — sin 77 cos 77 — 



(1 — cos 77) 



■n-'^-n 



* 



*, 



-(Vfet.,,.),.-^^^- ,.f^-^\, , (G23) 
V / 1 — C0S7; ^"'(l — cosTy)^ 

/— 3 sin 77 — 77 cos 77 — 277 2 + cos 7; 
yklM -r. r^ h -^^A/T- 7^, 



(1 — cost;)^ (1 — cos 77)2' 



.. = (Vfc/7 



A/ 



4 cos 

K/a/- 



+ {Dm) , 



3 sin 77 — 7/ cos Tj — 2ri 

(1 — COS7/)2 

77 + 77 sin 77(5 + COS 77) — 4 
(1 — cost;)^ 

Vkhiir]- sinr]) - A 

2 + cos rj 
(1 — cost;)^ 



(G24) 



- 4 sin^ 



7? 



sin 77(5 + cos 77) r 

-Dm— p. r^ Im{v- 

(1 — cos 771* 



'm 



. , Dm 
' sni77) — 



^/fc 



*, 



fc^/M 4 cos 7^ + 77 sin 77(5 - 



cos 7;) — 4 — 4 sin^ 

\ A 



(G25) 
77 



M (1- cos 77)4 

U^/'^Dm sin 77(5 + cos 77) 



*, 



M (1- cos 77)4 



(G26) 



M2(l-COS7/)5 



( — 2377— 1977 COS 77 



+277 sin 77 + 33 sin 77 + 9 sin 77 cos 77) 
PDm -19cos77 + 2sin277-23 



NP 



(1 — cos77)s 



(G27) 



From here on, the intermediate expressions become so 
large that we cannot reproduce them here; we only de- 
scribe how the calculation is done. We substitute (jG2ip 
- ((G27)) in (iGlSl) and multiply the resuh by (1 - cos 77) ^ 
to get a polynomial in (1 — cos 77) (of 6th degree). It is 
also a polynomial of 4th degree in 77. The function 77 is 
the only one in this polynomial that depends on t, the 
coefRcients of 77, (1 — cos 77) and their powers depend only 
on r. Thus we treat 77 and (1 — cos 77) as independent vari- 
ables, and the coefficients of their different powers must 
vanish separately. 

In W3 transformed in this way we now take the co- 
efficient of 77 , the resulting equation is a polynomial of 
degree 4 in (1 — cos 77): 



API 



M 



{e - kf 



[9{e /k -l)[l- cos rjf 

-78(£/fc-l)(l-cos77)3 

+ {2lQe/k - 276)(1 - 00877)^ 

-(189e/A: - 381)(1 - COS77) - 144] = 0. (G28) 



The term independent of (1 — cos 7;) clearly shows that 
Im = is a unique solution of this, independently of the 
value of k. 



So, we substitute Im = in the main large polynomial 
and in the resulting smaller polynomial we take the term 
independent of (1 — cos 77). The equation that results is: 



lAANpDM^/[k^{e-kf]^{). 



(G29) 



Here the unique solution is Dm = 0. But with /jv/ = 
Dm = we get * = from (|G24p . i.e. the Friedmann 
limit. Thus, 04^1 — ^4 = 0, which is one of the necessary 
conditions for the existence of RLPs, can in this case 
be fulfilled only when the Szekeres model trivializes to 
Friedmann. 

The calculation above was done for fc > 0. The calcu- 
lation with A: < is essentially the same and need not 
be done separately - it is enough to replace (fc, 77) in all 
equations with (— fc,i^. 

When fc = 0, we have necessarily e = +1 and the 
calculation must be done separately. Then we have: 



-(^ 



1/3 



{t-tsf^: 
1/3 



2 /9MV . x-4/3 



(G30) 



With the r-coordinate chosen so that M = M^r^ , where 
Mq is a constant, this simplifies W^ in (JGISI) to 



Wa 



dcf 



(64/9)MoViB/(t-iB)^ 



-(14/3) X iQ^'^M^^'^HB/ {t - tBr^°'^ 
+ (256/3)MoV5ts,,3 (i _ i^)-3 

+14 X m^'^Mo^'^Hs/ {t ~ ts)-"^ 



- ll2Ah' r^B./it-tBY 



(G31) 



-3 X 36i/3Mo^/ViB^,2 (^ _ ^^)-4/3 

+3 X m^'^Mo^l^r^tB,rrtB,r {t - ts)"'*^^ = 0. 

Now the independent variables are t and r, and t ap- 
pears always in the combination {t — ts)- Thus different 
powers of {t—ts) are linearly independent, and their coef- 
ficients must vanish separately. Whichever term we take, 
except for the last two, the result is always the same: 



tB^r = 



(G32) 



(because Mq = is the vacuum, i.e. Schwarzschild, limit 
of the L-T model). This guarantees that all the terms 
in (|G3ip vanish. However, as seen from (jG30p . tB,r = 
means Vf = 0, i.e. zero shear (see (|6.4I) and (j2.9|) '). i.e. 
the Friedmann limit. Thus, there are no non-radial RLPs 
also when fc = 0, which completes the proof in the case 
di^O. 

We go back to (|G15|) . where we assumed di ^ and 
proceed from there on to consider the case di = 0. In- 
stead of (|G18P we now get: 



W. 



-6$$„*,tV*2^2$$„*,«/* 
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e — k 

-3-^+3-^ =0. 



$ * 



(G33) 



As before, we begin by considering the case fc > 0. We 
multiply Wc, by ^^$(1 — cos 77)^ and substitute for $ and 
vf from (|G19p - (|G27p . What results is a polynomial of 
degree 6 in (1 — cos 77) and of degree 3 in 77. Taking the 
coefficient of rf we obtain: 



We 



MIm sin ry 



-6(e — fc)(l — cos 77)" 



+45(e - A:)(l - cosTy)^ - 81(e - fc)(l - cos?/) 
+30fc(l - cos 77) - 36/c] = 0. (G34) 

Looking at the term independent of (1 — cos 77) we see 
that the unique solution of this is Im — 0. 

So we substitute Im = in the main polynomial, and 
in the resulting expression we take the term independent 
of (1 — cosr;). The resulting equation is: 



36sinr;MDA/ 



(G35) 



The unique solution of this is Dm — 0, which, together 
with Im — 0, leads back to the Friedmann limit. Thus, 
no RLPs exist in nontrivial Szekeres spacetimes in this 
case, either. 

The argument given before, that the result for fc < 
follows by the substitution (fc, 77) — > (— fc, i^, is still valid. 
So we now consider di = with fc = 0. We substitute 
fc = (and, as is necessary, e — +1) in (|G33j) . multiply 
the result by ^, substitute then for ^I^ and $ from (|G30I) . 
and obtain: 

+2xiQ^'Hl^^'\tB/{t-tB)-'"' 
- {56 / 3) Marh B ^r {t-t By ^ 
+36^^^Mo^/^rtB,rr {t - IbY^'^ = 0. (G36) 

The coefficients of independent powers of {t — Ib) have to 
vanish separately, as explained before. Whichever term 
we take, except for the last one, the result is always the 
same: 



tB.r = 



(G37) 



and this implies the Friedmann limit in the same way 
as explained after (IG32p . This also guarantees that the 
whole of (|G36)) is fulfilled. 

Thus, in every case considered, the assumption that 
non-radial RLPs could exist leads to either the Fried- 
mann limit or the Schwarzschild limit. The final conclu- 
sion is that the only RLPs in the G3/S2 models are the 
radial null geodesies. D 



Appendix H: A detailed description of the model 
presented in Sec. IVIIII 



The algorithm used in the calculations discussed in Sec. 
IVIIII consists of following steps: 

1. First we set the observer at Rq (the present-day 
areal distance) and consider sources which are, at 
the present instant, away from the observer by 1 
Gly (« 306.6 Mpc). 

2. To calculate the evolution of the model one needs 
to follow the following points: 

• The radial coordinate is chosen to be the areal 
radius at the present instant: f — $(io,7')- 
However, to simplify the notation we will omit 
the bar and denote the new radial coordinate 
by r. 

• The chosen asymptotic cosmic background is 
an open Friedman model, i.e. flm — 0.3 and 
A = 0. The background density is then given 

by 



3iJ 

Pb = n,„. X Per = 0.3 X -^(1 + zf, (HI) 

where the Hubble constant is Hq = 72 km s~^ 
Mpc~^ 

• The initial time to is calculated from the fol- 
lowing formula for the background Friedman 
Universe 



tiz) 



Ho\l + I)-^dz 



ynr,,ii + zf + ii-a,n)ii + zf' 



(H2) 



• The age of the universe is assumed to be ev- 
erywhere the same: is = 0. 

• The function M{r) follows from (|2.5p . where 
the present-day density is 

p{to,r) ^po[l + S-Sexp (-rVa^)] 

• Because of the assumed spherical symmetry 
e" = 1. 

• The function k{r) can be calculated from 

dUD. 

• Then the evolution of the model can be calcu- 
lated from eq. 



3. We then find a null geodesic that joins the observer 
and the source. The angle between the direction 
towards the source and the direction towards the 
origin, at the present instant, is denoted as 7. 
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4. The null geodesies are found in the following man- 
ner: 

• Because of spherical symmetry we may set one 
of the angular components of the null vector 
to zero. We set fc"^ — 0. 

• The second angular component, fc^, follows 
from i?^fc^ = J — const = i?osin7, where 
Rq — R{to, ro), i.e. at the observer's position. 
This relation is a consequence of (|4.10|) and 
(14. lip and was derived in [l^l (see equation 
(3.26) in 111). 

• The radial component is evaluated from (|4.8p 
with Ei2 = = Ei3 = E,r, and E = 

• The time component of the null vector is found 
from k°'ka = 0. 



5. We then find two other null geodesies: one that 
will reach the observer in 1 Gy time, the other one 
that arrived at the observer's position 1 Gy ago. 
Because of the non-RLP effect these geodesies ap- 
proach the observer at angles that are different from 
7- 

6. The difference between these angles allows us to 
evaluate the rate of change of the angle 7. 
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